Abstract
Despite a characteristic indolent course, a substantial subset of follicular lymphoma (FL) patients has an early relapse with a poor outcome. Thus far, efforts to identify factors that predict survival have been unsuccessful. However, we and others have demonstrated the prognostic relevance of CD4+ T cells in the tumor microenvironment (TME) of FL and developed a bio-clinical risk model (called BioFLIPI) that offers improved risk stratification (Mondello et al, BCJ 2021). However, the mechanisms defining different BioFLIPI scores remain to be fully understood.
To identify the programs associated with the BioFLIPI scores, we first interrogated the RNA-seq of purified B cells from 67 newly diagnosed FL1-3A (BioFLIPI 1 n=21; BioFLIPI 2-3 n=36; BioFLIPI 4 n=10). We found 1,156 upregulated genes and downregulation of 930 genes in BioFLIPI 4 vs 1, including well-known immune molecules such as CD70, IL10RA, CD69, and CCR6 (FDR<0.05). Gene set enrichment analysis (GSEA) using the MSigDB confirmed a negative enrichment in multiple immune response and interferon gamma signatures and revealed striking positive enrichment for cell proliferation and antigen presentation pathways in BioFLIPI 4 compared to BioFLIPI 1 tumors. The same pathways were dysregulated in BioFLIPI 2-3 but a lower magnitude. To investigate the somatic mutations associated with different BioFLIPI scores, we performed WES of the same 67 FLs and further expanded this cohort for a total of 90 patients (BioFLIPI 1 n=26; BioFLIPI 2-3 n=51; BioFLIPI 4 n=13). We found enrichment for mutations in several regulators of the germinal center in BioFLIPI 4 and 2-3, such as MEF2B (p=0.029), FOXO1 (p=0.038) and IRF8 (p=0.035), suggesting the cooperation of multiple oncogenic events promoting lymphomagenesis. In contrast, BioFLIPI 1 tumors harbored a higher rate of mutations in TNFRSF14 gene, known to enhance anti-lymphoma immune response.
To investigate whether BioFLIPI modulates immune signatures of the microenvironment, we performed TME deconvolution by applying CIBERSORTx to bulk RNA-seq of these 67 FLs. We found that FLs with BioFLIPI 4 and 2-3 were enriched for CD4+ Tcm exh (p=0.02), Tem exh (p=0.04), Tfhexh (p=0.07) and NK cells (p=0.01) compared to BioFLIPI 1. In line with this data, CyTOF performed on a subset of these patients (BioFLIPI 1 n=12; BioFLIPI 2-3 n=25; BioFLIPI 4 n=6) confirmed progressive increase of CD4+ Tcm exh (p=0.006), Tem exh (p=0.009) and Tfh exh (p=0.025) cells from BioFLIPI 1 to BioFLIPI 4.
Single cell RNA/BCR-seq of 21 FLs (BioFLIPI 1 n=6; BioFLIPI 2-3 n=12; BioFLIPI 4 n=3) corroborated a higher abundance of Tem exh, cytotoxic (CTL) Tem exh, Tfhexh and NK cells in BioFLIPI 4 and 2-3. Conversely, CD8+ naïve T cells were prevalent in BioFLIPI 1 tumors. Differential gene expression in T cell clusters from BioFLIPI 4 vs BioFLIPI 1 FL revealed a marked upregulation of multiple immune genes involved in antigen presentation (e.g., HLA-DRB, HLA-DQA, HLA-DPA/B) and immune response (e.g., ETS1, TNFSF8, IL7R). GSEA showed a positive enrichment for antigen presentation signatures and negative enrichment for TNFα signaling and Oxphos metabolism. Similar findings were observed in BioFLIPI 2-3 FLs. Using CellChat analysis, we found a stronger interaction between malignant B cells and Tem exh and CTL Tem exh cells in BioFLIPI 4 and 2-3 vs BioFLIPI 1 tumors. This was linked to an increased strenght of MHCI-CD8 interaction and concordant decrease of MHCII-CD4. Notably, the suppressive BTLA-TNFRSF14 signaling axis was present in BioFLIPI 4 and 2-3 but not in BioFLIPI 1 FLs. This was accompanied by a decreased signaling through TNF/LTA-TNFRSF1B, CD80/CD86-CD28/CTLA4. Linking the BioFLIPI-specific immunogenicity to lymphomagenesis, we found that FL patients with BioFLIPI 4 had two-three folds higher risk of transformation to an aggressive lymphoma (12% vs 7%) or death (21% vs 3%) than those with BioFLIPI 1, suggesting that this immune dysfunction favors the expansion of a dominant tumor clone and promotes a more aggressive form of FL.
Collectively, our findings demonstrate that BioFLIPI scores in FL reflect distinct patterns of tumor immunogenicity that shape the TME. Notably, high BioFLIPI scores are paradoxically associated with enhanced antigen presentation alongside a suppressed immune response. This suggests a sustained antigen-specific T cell response that may drive immune exhaustion, promote tumor progression, and ultimately lead to histologic transformation